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1. The AI = 1 /2 rule 

One of the most striking hierarchies in hadronic physics is the famous AI = 1/2 rule, which 
refers to the experimental observation that the kaon decay amplitude in two pions with total isospin, 
/ = 0, is twenty times larger than that into ml = 2 state: 

T(K°^KK\ I=a ) =A a e iSa ^ =22.1. (1.1) 

The explanation of this large number is one of the notorious failures of large Nc which pre- 
dicts §\§ the ratio to be y/2. 

It was soon realized that there are many scales relevant in the dynamics of these decays, 
M\v,m c ,MK, ••• and maybe subleading orders in Nc could get enhanced by large renormalization 
group logarithms [ph. 

Below the scale of the W mass, this boson can be integrated out of the theory. The resulting 
effective Hamiltonian for the AS = 1 transitions can be obtained through the Operator Product 
Expansion (OPE) of two left-handed currents. The use of CPS symmetry and the flavour symmetry 
restricts the possible dimension six four-quark operators to just two, Qf and Qf, which are singlets 
under the SU (4)r and transform in the 84 and 20 representations of SU (4)i 

JifQ CD = 2VlG F {V us yV ud £ *f fif + k°Q°, (1.2) 



<T=± 



where 



fif = { ( s Yn p - u ) (uy^P-d) ± (sy^P-d) (uy^P-u) } - (u-> c), 
Qf=(m 2 u -mt){m d (sP + d)+m s {sP-d)y (1.3) 



The operators Qf do not contribute to the physical amplitudes and vanish identically if m u = m c . 
The 84 operator, Q\, contributes both to the AI = 1/2,3/2 transitions, while the 20, Q[, only 
contributes to AI = 1/2. This Hamiltonian can be run down to lower scales resumming leading 
logarithms through the renormalization group. A AI = 1/2 enhancement of the Wilson coefficients 
is then observed at low scales in commonly used schemes For example in the renormalization 
group invariant scheme (RGI) at 2-loops: 



*i (AO 



~2. (1.4) 

H=m c 



Although this goes in the right direction, a much bigger hierarchy must come from the hadronic 



matrix elements to match the experiment. Besides at /I = m c higher-order corrections to eq. (|L4j) 
are ^(100%). 

One can nevertheless try to push this perturbative analysis below the charm quark mass to 
obtain some qualitative understanding. When the charm is integrated out of the theory, one moves 
from a four-flavour theory to a three-flavour one. The OPE then allows any four-quark operator 
that transforms as the 27 or 8 dimensional representations of SU (3)l- Among the latter there are 
operators that are no longer the product of two left currents, but involve the product of left and right 
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currents or densities, such as the famous penguin operators. Arguments were put forward to argue 
that the hadronic matrix elements of the penguin operators are larger than those of the left current 
operators and this could be the origin of the rule [Qj. 

The computation of these matrix elements requires however a non-perturbative method. Many 
approaches to estimate the matrix elements using large Nc arguments or models have been pursued 
in the past. For a review and further references see e.g. Ref. []|]. Although plausible arguments 
seem to indicate that the rule could be the result of the accumulation of several instances of the 
octect enhacement, combined with large penguin matrix elements [Q], the approach relies on the 
use of perturbation theory down to dangerously low scales. 



2. The AI = 1/2 rule on the lattice 

It was realized in the early days of Lattice QCD [Q] that the AI = 1/2 rule would be a very 
well suited problem for the lattice approach, since it is such a large effect! Even if there are 
approximations like quenching, or in the presence of large systematic uncertainties, such a large 
enhancement would be hard to miss. 

In the pioneering work of Ref. m\, it was proposed to use the lattice to perform the matching 



of the AS = 1 effective Hamiltonian of eq. ( |1 .2| ) to an effective Hamiltonian in terms of the hadronic 
degrees of freedom, that is a chiral Lagrangian. The possibility to include the AS = 1 interactions in 
the chiral Lagrangian as a perturbation was first proposed by Cronin in Ref. In addition to the 
chiral Lagrangian that describes the strong interactions, Jz^pr, one would have an effective weak 
Hamiltonian with the same flavour symmetries as that in eq. (1.2), J%$. The operators Qf can be 



decomposed into 27 and 8 of SU (3)l- Therefore all the operators that can be constructed with the 
building blocks of the chiral Lagrangian, U,d^,M, with the same transformation properties as Qf 
should be included. At the leading order (LO) in a momentum or mass expansion there are just 
three of them, one 27-plet and two octects. The weak Hamiltonian can be writen [^] : 

^=2V2G f K rf c|^ 27 ^27 + 2g 8 ^8 + 2g / 8 ^|+H.c. , (2.1) 

where gs,g2i and g' & are low-energy couplings that contain the non-perturbative dynamics that is 
not fixed by symmetry arguments. The operators read 

^7 = [Kt^=\(m^ u +\m suud ) , (2.2) 



~~r ( d^UU j d^t/ 1 " , (2.3) 

4 V / ur \ / vs 

^ = \ I WAkU , (2-4) 

k=u,d,s 

& t = ^L(e w/N <MU + U Jf M<e- w / N! ) , (2.5) 

2 V / ds 

where we have made use of Tr [dfJJU^] = to simplify the expressions. 

If these couplings were known, the ratio of the AS = 1 amplitudes at LO in the chiral expansion 
would be: 

(2.6, 
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LO chiral perturbation theory is not expected to be very precise at the scale of Mk, but again the 
enhancement is such a large effect that as long as the expansion is reasonably well-behaved, the 
effect should already be there at the LO. 

It was then proposed in [j|] that these low-energy couplings could be determined by matching 
appropriate (the simplest) correlation functions between the chiral effective theory and lattice QCD, 
and in particular that this could be done through the computation of three-point functions and two- 
point functions. Note that to compute directly the amplitude K — > %%, four-point functions would 
be needed. 

This nice proposal turned out to be extremely difficult to implement in practice. Firstly the 
renormalization of four-fermion operators is extremely challenging when there is explicit breaking 
of chiral symmetry in the regularization. Not only there are a large number of additional mixings 
with wrong-chirality operators, but there is even mixing with lower dimensional ones and therefore 
power-diverging coefficients [10]. With the advent of Ginsparg -Wilson regularizations [ fTl] , |i2[ |l3| ], 



that preserve an exact chiral symmetry Q14Q, the renormalization of the effective Hamiltonian of 

ii. 



eq. (L2) becomes as simple as it is in the continuum JT5I] 1 



2± _ 7 ± .^±,bare . 7 ± ^)±,bare 
I — -1-^12^2 > 

Qf=Zf 2 Q^ e . (2.7) 



Indeed the computations of the AI = 1/2 rule in Refs. [|T(],[T7|] have been performed in the quenched 
approximation using domain-wall fermions that approximately satisfy the Ginsparg-Wilson 
relation. The lattice spacing used in these computations was however too low to keep the charm 
active, so they considered the effective Hamiltonian where the charm is integrated out perturba- 
tively. The effect of integrating out the charm brings in important complications. On the one hand, 
the renormalization involves power-divergent subtractions, which require a very good control over 



statistical as well as systematic uncertainties. A second difficulty was pointed out in Ref. Q19| ] and 
is related to the quenched ambiguities, which occur at the level of the OPE. More concretely there 
are operators in the OPE that are octects under the valence group and contain ghost quarks, in the 
supersymmetric formulation of the quenched approximation [EJj. Now, if these operators are not 
included in the OPE, by assuming for instance that the quenched approximation is only used to 
define the matrix elements at some low scale, it is unclear what the meaning of non-perturbative 
renormalization is, since it requires to compute the matrix elements in the quenched approximation 
up to very high scales. In other words it is not clear whether the use of the quenched approximation 
to compute matrix elements combined with the full theory OPE is really a consistent method. 

Finally the simulations in Refs. [T(| 17] were carried out at relatively large quark masses above 
m s /2. Chiral corrections were shown to be very large and a large systematic error resulted from 
the long chiral extrapolations. In the end, both collaborations found a large enhancement but there 
was almost a factor of 2 discrepancy between the two computations. 

3. New strategy 



The approach to the AI = 1/2 rule that was presented in Ref. [ 21 ] was designed not to reach the 



final result directly, but to try to reveal if the large enhancement is coming from one single leading 
1 Only in Ref. plf] the mixing with GIM suppressed operators has been discussed. 
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effect. The point is that on the lattice, the different physical scales that are involved in these decays 
can be modified at will in order to understand their relevance. In particular this is quite clearly the 
case for the charm quark mass. 

Most of the ideas that have been put forward to explain the enhancement are related in one 
way or the other to the charm quark mass. If the large enhancement is due to the large separation 
between m c and Aqcd or the up quark mass, there should be no enhancement whatsoever in a 
theory with a light charm quark. 

More concretely the strategy that we proposed to quantify the role of the charm quark mass is 
the following: 

• Step 1: light charm quark 

A theory with four degenerate quarks, m c = m u = = m s (GIM limit) is matched to a 
SU (4) chiral effective theory to extract the low-energy couplings that would mediate AS = 1 
transitions 

• Step 2: A x pj- 3> m c >»!„ = = m s 

If the charm gets significantly heavier than the up quark, but remains in the realm of the 
chiral theory, that is charmed mesons are well below the chiral theory cutoff 4nF , the charm 
quark can be integrated out from the SU (4) chiral theory to obtain the SU (3) chiral effective 
theory. If the charm quark is not too large this can be done analytically in chiral perturbation 
theory [^]. 

• Step 3: m c > A x p? 2> m u = ntd = m s 

If the charm gets too heavy to be describable in terms of an effective chiral theory, the match- 
ing to the SU (3) chiral theory has to be done non-perturbatively. The couplings g2i{m c ) and 
g%(m c ) can then be monitored as a function of the charm quark mass. 



4. K — > kk amplitudes in the GIM limit 

We will now describe the formulation of the problem in the GIM limit, that is for m u = rtid 



m s = m c . 



4.1 Lattice formulation 



In this limit there are only two operators in the OPE at first order in Gp , Qf in eq. (1.2), which 



renormalize therefore multiplicatively 



Qf 



y± ^±,bare 



(4.1) 



In the quenched approximation no spurious operator can appear [21]. 



It can be shown [21 ] that using overlap fermions [13, 14, 23] this simple renormalization pat 



tern is preserved. In Ref . [ 24 ] , the renormalization constants have been computed non-perturbatively 
through an intermediate matching to twisted-mass Wilson fermions at some large reference quark 
mass. The corresponding renormalization constants for twisted-mass Wilson fermions have been 
previously computed using the Schrodinger functional approach 



5 



The AI = 1/2 rule 



P. Hernandez 





bare P.T. 


MFI P.T. 


N.P 


Z\ i /z A 

^1 1 /^A 


0.525 
1.242 
0.657 


0.582 
1.193 
0.705 


0.584(62) 
1.15(12) 
0.561(61) 



Table 1: Comparison of the perturbative, mean-field improved and non-perturbative renormalization con- 
stants in the RGI scheme. 



The result for the renormalization constants in the RGI scheme (for details see [gjj, g4J) is 
summarized in Table [I], where the perturbative (one-loop) and mean-field estimates are also shown 
for comparison. 

4.2 %PT formulation 

In the GIM limit, the chiral Lagrangian has a SU(4) L x SU(4) R symetry group. It can be 
shown that at leading order in the momentum expansion, and in contrast with the SU (3) case, only 
two operators appear one 84 and a 20: 

je w xvr = 2V2G F (VusYVud £ g c [0 a ] 



where 



F 4 



4 

The normalization is such that in the large Nc limit 

[8 + }n c = [g-] Nc = I- (4-2) 
The ratio of the K —> Tin amplitudes in the GIM limit would be given by 

*>. = ±-(l + l£-) (43) 

and therefore at LO in the chiral expansion the hierarchy of the amplitudes is directly related to a 
hierarchy in the couplings g + and g_. Our primary goal is therefore to determine these couplings. 

We also note that no further operators appear in the quenched or partially quenched approxi- 
mation at this order, so there are no Golterman-Pallante ambiguities in the effective theory pi]]. 

4.3 The matching 

In order to determine the couplings g ± , we will match suitable correlation functions computed 
in lattice QCD to those computed to next-to-leading (NLO) order in the chiral expansion. Even 
though we will do the matching at very small quark masses, we will be in a finite volume, entering 
the £-regime of chiral perturbation theory p6[]. 
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We will consider three-point bare correlation functions of the four-fermion operators and two 
left-handed currents [27]: 

C a (x ,yo) = E<[/bWk» G ± (0) [Jw(y)}us) ■ (4.4) 

x,y 

More concretely we will match the ratio of these three-point functions to left-current correlators: 

R°{*o,yo) - jPSpl (4-5) 
C(jc )C(jo) 

where 

C(xo) =£<[/ioW]„p[/io(0)]p a > . (4.6) 

X 

Factorizing out the unknown couplings, the same ratio can be computed in the effective theory 

W a (xQ,y Q ) = , (4.7) 



n*oMyo) 



where 



^(xojo) = fd 3 x f d 3 y{[J w (x)) dll ^(0) [fw(y)]us) (4.8) 

^(xo) = | d 3 x([J LQ (x)] a pULo(0)]pa)- (4-9) 
The couplings can the be extracted from the matching: 



k°{M w ) n T " 7 \ , (4-10) 



where the Wilson coefficients are obtained in the same renormalization scheme as the operators, 
i.e. the RGI one. For the explicit expressions see [21]. 

4.4 @ c to NLO in %PT 

In order to obtain the LO couplings, the matching between lattice QCD and #PT should be 
done as close as possible to the chiral limit. The use of Ginsparg-Wilson regularization makes it 
possible to do simulations with extremely small quark masses: below a few MeV with volumes 
larger than 2 fm. However it becomes very costly to increase the volume much further than 2 fm 
or so. While the usual way to take the chiral limit is to first take V to infinity and only then take 
m to zero, there are advantages in taking m as small as possible at finite V, entering the so-called 
£ -regime defined by the condition: 

mLV<l. (4.11) 

For FL 3> 1, this limit is equivalent to ML <C 1, so the Compton wavelength of the pion is larger 
than the box size. Finite volume effects are large in this regime, however they are calculable in 
chiral perturbation theory [26]. The counting rules of the chiral expansion that ensure the condition 
of eq. ( 4.1 1 ) are: 

m£~ £ 4 LT\T- l ~s p~E (4.12) 
which are different to the usual ones. As a result of the new power-counting: 
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^-regime £-regime 



^qcd L 4 (D^D»U)(y) x 
L 5 (D^U' l D^U y) x 

u(y) 2 x 



Table 2: Strong and weak operators that would contribute to the observables considered in eqs.(4.9) at NLO 



in infinite-volume chiral perturbation theory. 5? = Ux^ +%U\ = i(Ux^ — XU r ), -S^ = Ud^U^md 



W^v = 2(c^Jz? v + d v y^), where % = 2mY./F 2 . The tensors tf- kj project onto the appropriate representations 



84 or 20. 



The zero-momentum modes of the pions become non-perturbative and have to be resumed to 
all orders. This is achieved by factorizing out the constant field configurations and treating 
them as collective variables in the partition functional. 

There is a reordering of the chiral expansion and at any given order less relevant couplings 
appear. In particular for the processes at hand the strong or weak higher-order couplings 
[30, 3l]] shown in Table ||] that would appear at NLO in infinite (or large enough) volume are 



all suppressed at the same order in the £ -regime. 

An additional complication is the quenched approximation, where we shall be working. The 
quenched chiral-perturbation-theory version of the £-regime was studied in detail in Ref. [[32|]. 
Analytical treatment is possible if averages are considered in fixed topological sectors. Ginsparg- 



Wilson fermions satisfy an exact index theorem [ 12] and therefore these averages can be computed 



on the lattice aswell. The different role of topology in the e and /^-regimes was first discused 



in [|29|]. The fixing of topology can be seen as a type of boundary condition, which should not 
affect the local properties of the theory, and therefore the dependence on the topological charge 
is in principle predictable in terms of the low-energy couplings of the effective theory, just as the 
finite-volume dependence is. 

The result of at NLO in the £-regime, in a fixed-topological sector of charge v, was 
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Figure 1: Left: NLO correction K as defined in eq. (4.13) as a function of the spatial extent of the box L 
in fm for three different ratios T/L. Right: 2R + (—T/3,T/3) for a lattice of 2 fm as a function of ML for 
T /L = 2. The band corresponds to varying the NLO couplings within a large reasonable range. The dotted 
and dashed lines are the °°-volume and the £-regime results respectively. 



computed in [28, 21] and the result is: 

l±K, (4.13) 



2 O*{xo,yo) = l ± [P~ 1/2 £i - Pk 

with p = T /L and j3i,&oo are shape coefficients of the box that depend only on p. The ratio turns 
out to be the same in the full and in the (partially-)quenched theories. It is quite remarkable that 
the result does not depend on the insertion of the currents XQ,yo, it is the same for all topological 
sectors v, and as expected no higher order coupling enters at this order. In the left plot of Figure |T] 
we show the numerical result for the NLO correction, K, as a function of the spatial extent of the 
box, for three values of the aspect ratio. If p is not too large 2 , the corrections are reasonably small 
for lattice sizes above 2 fm or so. 

Since we will simulate a number of quark masses, including those for which rriLV 3> 1, which 
is the usual regime of all previous calculations, we would need the results for the ratio & a in the 
/^-regime, where the counting rules are the same as in infinite volume and therefore at NLO the 
couplings of Table [2| become relevant. The result for 3% a in the /^-regime has been presented in 
[ 31 ] and is shown in the right plot of Figure |] as a function of ML for a lattice of extent 2 fm. 



too 



The band corresponds to changing the value of the unknown NLO couplings within a resonable 
range. For comparison the £-regime and oo-volume results are also shown. Surprisingly we find 
that deviations from the oo-volume expectations are significant up to ML < 5 for these observables. 

4.5 R a in lattice QCD 

The computation of R a for small quark masses in a finite volume is non-trivial. It was ob- 
served that large fluctuations in observables containing quark propagators occur in the kinematical 

2 The appropriate regime for large values of p is the so-called 5-regime, so actually the e-regime expansion breaks 
down in the p — > °o limit. 
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L/a 


T/a 


^low 


L[fm] 


HZ 


# cfgs 


e-regime 5.8485 


16 


32 


20 


2 


m. v /40,m,,/60 


0(800) 


jo-regime 5.8485 


16 


32 


20 


2 


m2/2 — m_ v /6 


0(200) 



Table 3: Simulation parameters 



conditions of the £-regime [33]. These fluctuations result from the fact that the low-lying spectrum 
of the Dirac operator is discrete in this regime. Both the low-lying eigenvalues of the massless 
Dirac operator and the splittings between them are controlled by the quantity (ZV) -1 : 

(A.) ~ ^ AA = -Xj ~ -^-i (4-14) 

therefore the eigenvalues of the massive Dirac operator in the £ -regime where m < (LV)^ 1 are of 
the same order as their splittings. In this situation, observables with point-to-all quark propagators 
can get large contributions from a few eigenfunctions. Space-time fluctuations in these eigenfunc- 
tions can lead to large and rare fluctuations in observables whenever peaks in the wavefunction 
happen to be near the fixed-point of the point-to-all propagator [fP|]. It has been shown [34] that 
these fluctuations could be cured provided rriLV ~ 1 through the technique of low-mode averaging 
(LMA) [34, 35]. The idea is to rewrite the quark propagator as a sum over its spectral decomposi- 
tion onto the subspace spanned by a few low-lying eigenvalues (low part) and the rest (heavy part). 
Once this representation is included in the correlation functions, those contributions that contain 
the low-parts only, can be averaged over all possible spatial insertions of the sources, because the 
wavefunctions of the low-modes are known in all points. The mixed contributions often can also 
be averaged provided a few additional inversions of the Dirac operator on the low -modes are per- 
formed. This is indeed the case for the two-point functions [ [341 ] and three-point functions [ 21 ] that 
we need in the GIM limit. The contributions from just the heavy parts remain unaveraged, but since 
the low-modes are no longer affecting these contributions, they should be much better behaved. 

The LMA typically requires the computation of the low-lying eigenvalues and eigenfunctions 
and the inversion of the Dirac operator on them, which is a non-negligible overhead, but still pays 
off as Figure ^ clearly shows. The two Montecarlo histories with and without LMA for one of the 
contractions of the three-point correlators of eq.(|4.4[) show that the improvement from LMA is very 
significant. 

In [36] we have presented the first results of a lattice determination of R a in the quenched 
approximation using overlap fermions [13]. The simulation parameters are summarized in Ta- 
ble |H The expected features of R° in the £-regime as predicted from chiral perturbation theory 
in eq. ( |4. 13 ) are well reproduced by the data. At large time separations, the ratio shows a flat be- 
haviour in xq and yo. There is no signal of v dependence in as shown in Figure ||. A weighted 
averaged is used to combine the result for all topologies. Finally there is no visible dependence 
with the quark mass for the small values considered. On the other hand the quark mass dependence 
of the ratios ^ in the /^-regime is quite significant. Note that in the /^-regime all topological sec- 
tors are averaged. In Figure || we show the results of the fits to NLO chiral perturbation theory 
expressions of the two combinations with smaller mass dependence: R + and the product R + R . 
The latter combination has the nice property that the NLO chiral corrections vanish at zero quark 
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MC history 



MC history (Q=3) 





Conf 



Figure 2: Left: Montecarlo history with and without LMA of the color-connected contraction in R a for the 
lightest mass in the ^-regime. Right: Montecarlo history with and without LMA of the color-disconnected 
contraction of R" in the e-regime for | v| =3. 
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Figure 3: Ratios as a function of the topological charge 2 < v < 10. 
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Figure 4: Left: NLO chiral perturbation theory fits of the bare R + as a funcion of the quark mass. Right: 
NLO chiral perturbation theory fits of the bare R + R as a function of the quark mass. The bands represent 
the statistical errors on the fitting function. 



mass. Each of the fits has two free parameters {{g + )bare, A + ) and {{g + g~)bare, A±), where A + ,A± 
is a combination of the NLO couplings that enter in the /^-regime. The bands contain the statistical 
(bootstrap) errors. Combining the results obtained from the fits of the bare ratios with the Wilson 
coefficients and the non-perturbative renormalization constants we obtain our main result: 

S+=0.51(3)(5)(6), g-=2.6(l)(3)(3), (4.15) 

where the first error is statistical, the second comes from the renormalization factor and the third is 
a systematic error estimated from the dispersion of the values of g ± obtained from different fiting 
strategies (for example fitting only the £-regime or only the ^-regime points). Although the results 
have been obtained for just one lattice spacing, scaling studies of several observables with overlap 



fermions have shown that scaling violations tend to be very small [ ]37[ ] . 

These numbers can be compared with those that would be obtained from experiment if the 
A/ = 3/2 and AI = 1/2 amplitudes would be matched to the corresponding ones in LO chiral 
perturbation theory in the GIM limit: 

g + ~0.5, g-~10.4. (4.16) 

Therefore the value of g + is strikingly close to experiment already in the GIM limit, but the value 
of g_ is a factor ~ 4 smaller. A significant enhancement is therefore observed in this limit 

^-6, (4.17) 

which is not large enough to explain the experimental ratio but is already significant and cannot be 
abscribed to penguin operators nor penguin contractions. 
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5. Towards a heavier charm 



Once the LECS in the GIM limit are known, the charm quark can be increased. If the charm 
is still light so that charmed mesons can be treated as chiral degrees of freedom, it is possible to 
match the two chiral theories with and without the charmed mesons analytically, or in other words 
the charmed mesons can be integrated out of the SU (4) chiral theory. In this way one recovers an 
SU (3) chiral Lagrangian where the couplings g2i and can be computed as functions of g± and 
the charm quark mass. Only when the charm quark mass satisfies 



m,X mJL 2 
< -j^ < (4nF) 2 . 



F 2 



is this perturbative matching reliable. 



This exercise was carried out in Refs. [22, 21] with the result: 



(5.1) 



gzj(m c ) = -g 



g + 1 + 15 



-^Lln^ 

(4k F) 2 M c 



+ g~ 1 + 3 



(4nF) 2 M c 



(5.2) 
(5.3) 



where A x contains the unknown information on the NLO couplings that are relevant as m c in- 
creases. In Figure § the dependence on m c of the ratio gz/gn is shown for two values of the 
ratio of couplings in the GIM limit g~ /g + = 1,5. The surprising observation is that only the oc- 
tect coupling has a logarithmic enhancement. Unfortunately NLO couplings are needed to have 
predictability. The effect of these is represented by the band where they have been varied within 
reasonable values (the associated scale A x = 1 — 4 GeV). For a ratio of g~ /g + ~ 5, close to the 
value we have obtained in the GIM limit, there could be a large effect coming from the integration 
of the charm. Unfortunately there is no much predictability unless the NLO couplings are known. 
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In order to go beyond chiral perturbation theory it is necessary to do a non-perturbative matching 
this time with a heavy charm quark, which is the next step of our project. 

It is well known that the case of a heavier charm will bring additional challenges. On the 
numerical side, the computation of three-point functions with a heavy charm requires the evaluation 
of the penguin contractions, which involve a point-point propagator. We are confident that LMA 
will also help in this case, but this has yet to be demonstrated. On the theoretical side, the quenched 
ambiguities of Golterman and Pallante will be present in the SU (3) chiral effective theory. It has 
been shown in pT| ] that these quenching ambiguities are rather mild in the £-regime, and it is in 
principle possible to disentangle the "physical" couplings from the spurious ones. 



6. Conclusions 

The AI = 1/2 rule remains a big challenge for lattice QCD. We have presented a well-defined 
strategy to quantify the role of the charm quark mass and in particular of the penguin contractions 



in the enhancement Q2 1[] . The idea is to compute the low-energy couplings mediating the AS = 1 
transitions in a theory where the charm is light and degenerate with the remaining three quarks 
(GIM limit) and compare them with those in a theory with just three light flavours, where the 
corresponding couplings can be computed as a function of the charm quark mass. 

The GIM limit is easier to treat in many respects and the first results for the couplings have re- 



cently been presented in Ref . [ 36 ] in the quenched approximation. The low-energy couplings in this 
limit already show a significant enhancement of the AI = 1/2 type. Even though the enhancement 
is not large enough to match the experiment, it already indicates that penguin operator/contractions 
cannot be the whole story. 
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